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Abstract 

We build a statistical ensemble representation of two econoniic 
models describing respectively, in simplified terms, a payment system 
and a credit market. To this purpose we adopt the Boltzmann-Gibbs 
distribution where the role of the Hamiltonian is taken by the total 
money supply (i.e. including money created from debt) of a set of 
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interacting economic agents. As a result, we can read the main ther- 
modynamic quantities in terms of monetary ones. In particular, we 
define for the credit market model a work term which is related to the 
impact of monetary policy on credit creation. Furthermore, with our 
formalism we recover and extend some results concerning the temper- 
ature of an economic system, previously presented in the literature by 
considering only the monetary base as conserved quantity. Finally, we 
study the statistical ensemble for the Pareto distribution. 

PACS Numbers: 89.65 Gh, 89.75 Da, 87.23 Ge 

Keywords: Money, Debt, Thermodynamics, Economic equilibrium. 

1 Introduction 

Historically, the first attempt to establish a link between classical thermody- 
namics and economics is due to Samuelson [I], who remained nevertheless 
quite skeptical about the possibility to introduce concepts borrowed from 
statistical mechanics, like entropy, into economic^. This is not surprising, 
since the neoclassical notion of economic equilibrium, in analogy with clas- 
sical mechanics, refers to a single point in the phase space, which is the 
solution of a maximization problem, and not to an equilibrium probability 
distribution like in the case of statistical mechanics. 

When heterodox economists criticize mainstream economic theory by un- 
derlining that the economy is not a system in equilibrium, they refer mainly 
to the first, static, notion of equilibrium. In fact, the general idea of employ- 
ing statistical equilibrium as a key tool for economics is advocated by some 
leading heterodox economists (see e. g. [6]), while models borrowed from 
statistical physics have been successfully applied to microeconomic models 
where agents are not perfect optimizers and/or not perfectly informed, like 
for instance in [7] or [8]. On the other hand, unconventional economists 
are rightly worried about the risk in substitute the unrealistic neoclassi- 
cal paradigm with an equally unrealistic paradigm derived from statistical 
physics [9]. In fact, it is quite unreasonable to think of long run economic 
evolution, which is the outcome of a complex out-of-equilibrium dynamics, 
as being explained by conservation laws, although we easily admit that such 
an economic world, described by stable probability distributions, would be 
indeed one in which our daily life would be much easier. 

mapping between classical thermodynamics and economics can be found 
in ^ , without reference to the underlying statistical structure (see also [3l |4l 

m)- 
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Some authors [lOj argue that if we accept the possibiUty that 'over a hm- 
ited period of time, an economic system might behave as though it were 
quasi equihbrate' in a statistical sense, the knowledge of equilibrium states 
for economic systems would then explain their possible 'thermodynamical' 
evolution. Unfortunately, from the perspective of economists this otherwise 
attractive bargain is based on a claim which is still overreaching, since a na- 
tional economy is still too complex to behave conservatively even in the very 
short run. We should be much more modest since, as explained by [9], the 
only conservative process in the economic sphere is exchange which, apart 
from being a short run phenomenon, delivers an inherently partial represen- 
tation of the economic reality. In particular, equilibrium distributions are 
conceivable only by separating temporally production and exchange, in such 
a way that supply is kept fixed over the time span considered in the model. 
Even taking into account this limitation, we can think of many economic 
contexts in which the application of the conservation principle could turn 
out to be useful. In general, conservation may be assumed to hold in all 
those markets where supply may be taken to change slowly, while exchange 
takes place at a relatively fast rate, enabling a relaxation of the system to- 
wards statistical equilibrium. These markets are quite common in the eco- 
nomic reality. A first simple example is provided by the daily functioning of 
fish markets, since fishermen cannot expand or shrink their supply of fresh 
fish once the market is open. A second, perhaps more relevant, example is 
provided by the real estate market, since the supply of new houses comes 
necessarily with a significant time lag with respect to positive variations of 
demand, and a negative adjustment of supply is generally impossible under 
normal conditions. A third example is provided by financial markets, where 
the supply of financial assets (stocks, bonds, etc.) can be easily assumed 
to be fixed on a daily basis. A fourth, less obvious, example is the credit 
market, since the supply of currency, customary referred to as the monetary 
base, is controlled by the central bank, and can be assumed to be fixed over 
a short time span, whereas money supply is determined endogenously by 
agents interactions through lending and borrowing. 

Recently some authors have followed a more ambitious path, by tackling a 
macroeconomic problem, namely the determination of the most likely data 
generating process for the empirical distribution of income and wealth at the 
national level [HI Hills], with the tools of equilibrium statistical mechanics. 
In these papers it is argued that the main ingredient for the derivation of the 
Boltzmann-Gibbs law is the existence of a conserved quantity [14J. In this 
context, the probability distribution for an energy value e is P(e) ~ e"*^/^, 
where T is the temperature. The authors of jll^ [T2| [T3] argue that money, 
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and not material wealth, could act as a conserved quantity in a closed eco- 
nomic system. In fact, material wealth as goods and properties can be 
consumed and destroyed, while an agent (also a firm) cannot print money, 
but only exchange it with other agents. 

This assumption has been criticized by other authors [TS], who have un- 
derlined that money is not a conserved quantity since, as underlined above, 
the endogenous lending process generally changes the amount of money in 
the system. From this perspective, the introduction of debt under the form 
of negative cash [11] clashes with economic realism since cash by definition 
cannot be negative. In order to be more realistic, we need our definition of 
money to be consistent with the well established standard notions of mon- 
etary base and money supply. In order to keep terminology as simple as 
possible, in this paper we choose to define money (cash), unless otherwise 
specified, as synonym of the monetary base, while the sum of currency and 
of highly liquid credit assets, like demand deposits (which we will indicate 
simply as 'credit'), will be labeled as total money supply. 

Furthermore, in [11] a Boltzmann-Gibbs equilibrium distribution for 
money is obtained as a result of numerical agent-based simulations. A sim- 
ilar study, with an attempt to introduce an ensemble for the conservation 
of money, is given in [16]. In fact, the introduction of a Boltzmann-Gibbs 
distribution, in order to build a statistical mechanics for money, should also 
be justified by introducing an ensemble starting from the microscopic eco- 
nomic variables from which total money supply is composed, similarly to 
the introduction of the coordinates for a given ensemble. In this way, the 
introduction of aggregate economic variables is obtained after an integration 
over the microscopic variables which describe the total money supply of a 
given economic system. 

This approach is the starting point of our paper. We will see that this 
representation adapts naturally to economic interpretation. In particular, 
with this formalism it is possible to introduce credit as a further variable 
with respect to the conserved monetary base, providing results which are 
consistent with the endogeneity of total money supply required for a sound 
economic interpretation. 

The plan of the paper is the following: in the next section we introduce 
our formalism together with its economic interpretation, in section 3,4,5 we 
introduce the statistical ensembles for money, in section 6 we analyze some 
macroeconomic relations suggested by the thermodynamics analogy while 
in section 7 we track some conclusions and final remarks. The appendix is 
devoted to the application of our method to the Pareto law. 
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2 Preliminaries 



In this paper we want to introduce sound statistical ensembles for econom- 
ical systems with money (cash) and credit-debt. First of all, in our model 
the particles of the standard physical approach are substituted by N 
interacting economical agents. In this context, we use the coordinate Xi 
to label the cash (money) possessed by the ith agent, while yi represents 
its credit-debt variable. In our context, we have no need to employ the 
complete Hamiltonian formalism of the ordinary statistical mechanics. In 
particular, the variables {xi,yi} we use are not conjugated in the sense of 
the Hamiltonian mechanics. To build an ensemble we have only need of a 
conservation law. We take as a conserved quantity the total money supply 
function M, which can be considered constant over a suitable time interval. 
Moreover, the function M will be a function of the chosen coordinates, i.e. 
M = M{xi,yi). 

At this point the ergodic hypothesis comes into action which permits us, 
given a function f{xi,yi), to express its average with respect to the time in 
terms of an average over the ensemble at fixed M: 

1= f{x,y)p{x,y)dxdy, (1) 

J M=const 

where p{x, y) denotes the probability distribution of the ensemble. Taking 
into account these considerations, our fundamental aim in this paper is to 
show that ensembles can be introduced in a consistent way for an economic 
system. To this purpose, we study simple money functions reproducing 
known results present in the literature, such as the following: 

N 

M = ^ (x, + y,) (2) 

i=l 

In expression ^ both Xi and yi can be expressed in monetary units. We 
stress that the expression ([2]) is the simplest choice we can consider. More 
sophisticated models can be introduced by means of an interacting term 
I{x, y) in the money supply function M. Once the ensembles have been de- 
fined, we can introduce in a consistent way all the thermodynamic functions 
in terms of and of economic variables, and finally we can study economic 
(thermodynamical) transformations . 
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3 Microcanonical ensemble 



In order to introduce the microcanonical ensemble we consider an isolated 
economic system with N agents, where money M is fixed. As in the usual 
microcanonical ensemble in statistical mechanics, if we integrate over all the 
available volume of the configuration space spanned by the variables {x, y} 
with M = m = constant (the overline denotes the average over the whole 
configuration space), then we have O = fj^^^d^x d^y = (see jl7] and 
references therein). As usual in statistical mechanics, we introduce a thick 
shell A where A << m and define 

r(»)^/ (3) 

Jni<M<m+A 

where k is a normalization factor such that T is dimensionless. In the fol- 
lowing we set k = 1. The average of a given quantity / is given by 

7^^. (4) 

To calculate r(m) it is convenient to introduce the following integral: 

S(m) = / d^x d^y, (5) 

JM<m 

where 

r(m) = S(m + A) - S(m) ~ --^A. (6) 

dm 

The functional r(m) measures the number of microscopic realizations of a 
given economic system. The entropy S can be defined in the usual way: 

5 = Inr(m) = InS(m), (7) 

where the last equality in ([7]) follows for > > 1 (see [T^ ) . In our descrip- 
tion of an economic system, entropy keeps its ordinary statistical meaning, 
i.e. it is proportional to the number of micro configurations of the system, 
which is consistent with a given distribution of the A^ agents with respect 
to X and y. Then, by analogy with ordinary statistical systems, the equilib- 
rium configuration is the one that maximizes the entropy. This condition, 
together with the conservation of money (energy) leads to the ensemble dis- 
tribution. 
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After writing S = S{m,V), where Vy = J dy, we can define the analogues 
in our context of temperature T and pressure P: 

Furthermore, we define the analogue of the free Helmholtz energy F: 

F = m-TS, (10) 

dF = dm- TdS - SdT. (11) 

From dS]) and ([9]) we get: 

TdS = dm + PdV. (12) 

Equation (jl2p is the analogue of the first thermodynamic principle. From 
PT]) and we have: 

dF = -PdV - SdT, (13) 

In the following section we give an economic interpretation of the quanti- 
ties P and T, while Vy is naturally interpreted as the maximal amount of 
credit/debt per agent allowed in the system. 

As a first warm-up example, we want to reproduce known result in the 
literature [IT], 

N 

M = ^x, (15) 
1=1 

In the expression (jl5p we have supposed one component for any pair of 
coordinates (x,y), with x G [0,oo). This is the simplest assumption we can 
make. Integration of ([5]) is easy and it gives: 

N 

= (16) 



By using ([7|) and ([9]) we have: 



777 777 A/^'T"' 

S = Nilnm + lnVy)-lnNl, T = - P = — = —. (17) 

iV Vy Vy 
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The expression for T of ()17p is the same as that of the model of Dragulescu 
and Yakovenko [11]. In the form ()17p . T measures the mean money per 
agent. In the foUowing section we obtain different models in terms of the 
more usual canonical ensemble. 

4 Canonical ensemble 

As the usual canonical ensemble, we have an economic system with A'^i 
agents with money mi interacting with another system with agents (with 
^^^^2 >> -^i) with money m2. Under the assumptions that total money 
m = mi + 1712 is conserved and that: 

r(m) ~ ri(mi)r2(m2 = m — mi), (18) 

we can calculate the probability distribution at equilibrium of the subsys- 
tem 1 independently of the behavior of reservoir 2. As a consequence, for 
an economic system with agents, the money function M and economic 
temperature T are given by 

dP=^^-^^^^^, (19) 
Z = J d'^xd^'ye-^, (20) 

Moreover, by a simple algebra we obtain 

Z = e'^, (21) 

with F given by (jlOp . It is natural to introduce the analogue of the Maxwell- 
Boltzmann distribution for a single agent. Actually the infinitesimal proba- 
bility dPi that the ith agent is in a point of its configurational space {x,y} 
independently of the behavior of the other agents can be written as 

dPi{xi,yi) = ce~~dxidyi, (22) 

where c is a normalization constant. Note that also in |1H \T2\ [T3] a dis- 
tribution similar to is considered, but integrated in a different space. 
In fact, in [TTl |12l [13] the exponential in ([22]) is integrated with respect to 
drrii, while in (|22p according to the definition of the Boltzmann-Maxwell 
distribution for a classical ideal gas, it is integrated with respect to {xi,yi}. 
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4.1 Model of a payment system with only cash 

We first derive the model (jl5|) in terms of the usual Lagrange multipliers 
(LM) by a discretization of the system and after we translate the model in 
terms of ensembles. The model describes the following economic situation. 
The generic agent j with initial cash Xj pays a certain amount of money 
Ax to another agent k with initial cash x^. In the final configuration we 
have the agent j with cash = Xj — Ax and the agent k with final cash 
xj/ = x^. + Ax with the sum conserved. This situation, as underlined by 
[18], describes a payment system, which is related to a unobservable flow of 
goods and services across agents. 

The number of configurations T of agents with money is 



r(K}) 



ni\n2\ • • • rifc! 



(23) 



The entropy S = lnr({nfc}) must be maximized with the LM a and /3 and 
with the constraints 

n{p) = N, ^ n{x)x = m= (M) , (24) 

p X 

by taking 

S = S -a^n{x) - ^^n{x)x. (25) 

X X 

The maximization of (j25p gives n{x) = e~"~^^. By introducing this ex- 
pression in (I24p and changing summation with integration, i.e. ^ n{x) — )• 
/ n{x)dx dy, we obtain 1//3 = T = m/N. 
Within the canonical ensemble we have 



AT 



d'^xd'^ye — —. (26) 

We easily obtain, by a factorization of (|26j) . Z = V^T^ . By means of 
(j2ip we have F = —TlnZ. By means of (jl4p and (jlOp we obtain for T 
the expression T = m/N. This case is similar to the perfect gas one of 
the usual statistical mechanics. The gas is in a box of volume Vy. As 
stated above, in the economic system, Vy can be interpreted as the maximal 
credit level allowed in the system, although in this case the function M does 
not contain explicitly credit-debt terms. Any thermodynamical economical 
transformation which increases Vy simply increases the allowed credit level. 
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4.2 Models of a payment system with debt 

In the model [ llj debt is introduced by allowing for negative money. In 
our approach money is always positive, and as a consequence our analogue 

N 

of [TT] is given by M = ^ y^, where Ui € [—d,oo),d > 0. The economic 

1=1 

dynamics of this model is the same as that of the previous subsection above, 
but y here represents a current account balance with maximal overdraft —d. 
Then, if 7/j > agent i has a credit claim towards a bank, supposed to be 
external to the system, which he can use to pay for goods and services, while 
if j/i < the same agent is borrowing from the same bank in order to fulfill 
the same payments. In terms of the LM technique we have Yl '^(y) = ^ 
and Yl''^{y)y = Qo = {^)- By the same technique as that of the previous 
section , with ^ n(y) J n{y)dy dx and y G [—d, oo) we obtain 

T=^ + d. (27) 
In terms of our formalism we have 

JV 

d'^xd^'ye-'^. (28) 

By performing the same calculations as that of the previous subsection we 
get Z = V^T^e~ . As a result we obtain for T the expression (1271) . In 
this case Vx represents the maximal amount of cash allowed in the system. 
Note that we can also have Qo < 0, with the obvious condition T > 0, i.e. 
d > —Qo/N. Moreover, we can build a model with d < and the formula 
(|23p is again valid. The model (j28p can be generalized further. In fact, 
suppose that each agent i has different current accounts {yii,yi2 ■ ■ -Uiri} 
with maximal overdrafts {dn, di2 ■ ■ ■ dir-}. Then we have 

^ = ^ + E = + (d) (29) 

where R = X^i^i^^i- We now consider the more general money function 
given by which represents a system where agents make payments either 
by cash or through a single current account (e.g. by using bank checks or 
credit cards). It should be noticed that in this case the constraint is on 
rrii = Xi + yi, since we suppose that agents can freely convert Xi into yi and 
vice versa. As a result, no constraints are imposed on X^yj or X^Xj. 
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In the first place, we can build a model where both credit and debt are 
allowed. In terms of LM we have ^ n{m) = N, n[m)m = (M) where 
Hi G [— d, (X)),(i > 0. By means of the same technique we previously used 
(remember that in this case ^ n(m) — )■ J n[ra)dx dy with (j23p again valid) 
we obtain 

T=5(^ + <i). (30) 

In order to translate the model in terms of the ensemble formalism, we 
proceed exactly as in subsections 4.1 and 4.2. As a result of the integration 
of d^ni) and thanks to ([TU]) and ([HD we obtain 

F = -2NT\n{T) - Nd (31) 
S = 2N\n{T) + 2N (32) 

T ^ i(^ + .). (33) 

We remark that the temperature we get is one half the one obtained in [11] , 
This factor is related to the two alternative payment means for each agent 
we have in our model. In fact, it's easy to see that (|33p is a special case of 
([29]) with yn = Xi, yi2 = yi, Vi = const = 2, d^ = Q and dy = const = d. 
Thus, the introduction of alternative payments means for a fixed m makes 
the tail of the Boltzmann distribution for money thinner, by decreasing the 
value of T. 

We now analyze a model, defined by the same money function ([2|) , where 
credit is not allowed, i.e. yi G [— 0] . This model is consistent, although not 
necessarily realistic, if we allow total cash to be increased by debt creation. 
To this purpose, consider an initial condition in which a generic agent '1' has 
x\ = xi^ 2/1=0 and another labeled '2' has X2 = X2, y2= 0- In a first phase 
T the agent '2' borrows money from '1'. We have x\ = yi — Ax, y[ = Ax = 
Ay and X2 = X2 + Ax, = Moreover, the agent '1' can, by means of 

a central bank, change credit with respect to the agent '2' with money cash. 
Hence, in the final state 'IT we have x{^ = xi — Ax + Ax = xi, y\^ = and 
X2^ = X2, 2/2^ = 2/2- ^^^^ way, the total amount of cash is increased, since 
xl^ + = xi + X2 + Ax, with a net creation of debt —Ax = Ay. This 
result can be interpreted as a debt of the agent with a central bank. This 
model can be solved by integrating with yi £ [—d,0]. 

By integrating equation (j26p for Z and with the help of (j2ip . (jlOp and (jl4p 
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we get: 

F = -2iVrin(r) - iVT In (e*' - 1) , (34) 



S = 2Nln{T) + 2N + N\n(eT -1) " (35) 

^ ^ T eT-V 



^ = ^NT-— ^. (36) 

eT — I 



Note that the expression m in (j36p is much more involved than the one of 
However, for T » d expression ()36p for the temperature T becomes 

r=^ + ^+o(d), (37) 

which is a result which differs, at the given approximation order, from that 
in only in the coefficient 1/2 (instead of 1) of the linear term d. 



4.3 Model of a credit market 

A credit market differs from a payment system in that it is constrained 
to respect the basic accounting identity over interactions between pairs of 
agents. If the accounting of a credit relationship is correctly reproduced, 
when agent 1 lends to agent 2 we have 

Ml = (xi - Am) + (di + Am) = consti (38) 
M2 = (x2 + Am) + (d2 - Am) = const2 (39) 

In this case, di stands for net debt, i.e. di = assets^ — liabilities j. In our case, 
we have a single class of assets (credit) and of liabilities (debt) respectively. 
Furthermore, the following additional constraints must hold: 

T.X^ = MQ (40) 

^ assetsj — Yli liabilitiesj = Qq (41) 

The first constraint means that the monetary base Mq is fixed, i.e. that 
agents, by exchanging credit, can increase or decrease only the overall money 
supply (i. e. the sum of money (cash) and credit), while Mq remains fixed. 
The second constraint means that the system is closed, i.e. any credit claim 
of some agent i in the system represents a debt for some other agent j in 
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the system and vice versa. In the fohowing, we suppose without loss of 
generahty that Qq = 0, so that credit and debt distributions are identical 
The money function for this model is given by ^ = m, where yi stands 
for credit, i.e. yi G [0,+oo]. In fact, since Mq is fixed, the only source of 
fluctuations for M in the canonical ensemble is given by credit /debt y,. 
This fact allows us to obtain a well defined temperature for the system. In 
terms of the ensemble formalism we have the following: 

N 

E Vi 



ry e-^ = j d'^'y e"— = M^'r\ (42) 

where in (j42p we have posed Mq = Vx- This assumption can be understood 
in the following way. Since Mq is fixed, we have that X] ~ -^O) and 
changing summation with integration we have Jj2x <Mo'^'''" ~ = Mq. 
Furthermore, it is easy to see that the expression (I17p for P still holds, with 
V = Vx = Mq. By following the same steps as the previous section we obtain 
again 

T=^, (43) 

Of course, we could obtain the same expression by integrating only over the 
yi G [0, +oo]. Nevertheless, we maintain the usual representation in terms of 
Xi and yi in order to reproduce correctly the main thermodynamical equa- 
tions. In fact, as we will see, the variations of ^ = Mq as determined 
by monetary policy are needed in order to introduce a "work" term in the 
variation of credit. In this context, the variations of m reflect instead a 
change of attitude on the part of the agents, which may become more or less 
oriented towards credit creation. 

Finally, we introduce a generalization of the money function ([2]) of the fol- 
lowing kind: 

TV I 

^ = J2J2y^^ (44) 
i=i j=i 

where the index i runs over the agents while j runs over the different asset 
classes allowed in the model. This generalization turns out to be useful to 
describe a full fledged flnancial market with different kinds of assets and 
liabilities. In this case we obtain the following 

F = -TNIlnT (45) 
S = NIlnT + NI (46) 



13 



Thus, if one allows more kinds of assets/liabilities into the system, the likeli- 
hood of observing large values of rrii, as well as large deviations from m, will 
be lower, since T is lower for the same m with respect to the one-dimensional 
case with temperature (03]). In other words, diversification makes the sys- 
tem more stable, a result which is consistent with standard economic theory. 
For a further discussion on debt the reader can refer also to [18\ [T9] . In par- 
ticular, in the model present in [19] a constraint on debt is introduced by 
adopting the fractional reserve system of banking. Even this model could 
be reproduced with our formalism in terms of a statistical ensemble. 

5 Grand Canonical Ensemble 

This kind of ensemble is aimed at describing the probability distribution 
of a system in which the agent number is not conserved. The probability 
distribution with N agents is given by 

dP^j, = e(-^+^-f d^y, (48) 
where is the analogue of the chemical potential and: 

dS=^ + ^dV- ^dN, (49) 

dF{T, V, N) = -PdV - SdT + fidN, (50) 
dF\ fdm\ , , 

The last equality in (|5ip follows from the Maxwell relations [H] . From (|49p 
we have: 

TdS = dm + PdV - ^dN. (52) 
The analogue of the Gibbs-Duhem relation gives 

SdT = VdP - Ndfi. (53) 

In this context, the 'potential' term /x takes into account the contribution to 
m caused by a variation of the number of the agents. For an interpretation 
of fi in terms of immigration see |20j . 



6 Thermodynamics of Money: Some Relations 

The starting point for a thermodynamics of money is provided by equa- 
tion ()52p . In the following, we will refer to the credit market model as a 
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benchmark to provide a sound economic interpretation of thermodynamical 
relationships. We should remember that in this model V = Vx, i.e. V is 
related to the cash variable x and thus to the monetary base. First of all, 
thanks to ([5T|) and by taking (5, V) = constant we get 



Hence, the analogue of the thermodynamic chemical potential in our context 
has to be identified with the variation of credit induced by a variation of the 
agent number. In complete analogy with thermodynamics we set TdS = dC 
in ()52p . where dC represents the increase (or decrease) of credit, which is 
caused respectively by a change in the attitude of agents (dm), by a change 
of the monetary base induced by the central hsi\k[PdV), or by a change 
in the number of agents {^dN). In fact, further pursuing the analogy with 
thermodynamics we can compute the quantity 



L,^f = / PdV, (55) 

where i, / refer to an initial and a final economic state characterized by 
{T, P,V, N). The integral (I55p depends in general on the path from i to 
/. Moreover, as underlined, V = J dx is a measure of the monetary base 
present in the system. Then a change in V represents a change of monetary 
policy by the central bank, and it can be naturally interpreted as the work 
done on the system by the central bank itself. Summarizing, for quasi-static 
economic transformations, we can write 

dC = dm + dL- fidN, (56) 

In this way, equation (I54p can be rewritten as: 

Furthermore, if the economic transformations can be considered as quasi- 
static, then we can insert in (j55p the expression (jl7p for P which is valid 
for a 'free' agent. In our case, P = y^, i.e. P in our context captures 
the relationship between credit and money, by describing the magnitude of 
credit for a given value of the monetary base. We remark that, if we followed 

the fractional reserve system model of [19], we should write m = ^ 1^ V, 



15 



where r is again controlled by the central bank. In this case, credit would 
be completely determined by the central bank itself. 

If we consider more general transformations than quasi-static ones, equation 
([56]) becomes: 

6C = dm + 6L- fidN, (58) 

where the non exact differential 6 takes the place of the exact one d. In this 
way we can consider economic transformation that are not 'quasi-static'. 
As an example we can study the analogue of an irreversible transformation 
in the usual thermodynamics by considering an isolated economic system 
evolving 'spontaneously'. In this case, as the system is isolated, {C,L,N) 
are left unchanged. Hence also dm = 0, and the economic temperature 
remains unchanged (as it happens for a free expansion of an isolated gas), 
provided that the money function M contains only 'free' terms yi. We can 
also calculate the entropy S as: 

where the integral (j59|) must be calculated over any quasi-static transfor- 
mation connecting the initial and the final economic states. Note that the 
formula (|58p . for a change with fixed C and agent number N, yields: 

dm + 5L = 0, (60) 

Adiabatic transformations are then naturally interpreted as those where the 
central bank counteracts the endogenous increase (decrease) of credit with 
an opposite variation of the monetary base in such a way that the overall 
credit remains fixed. 

We have seen that in the credit market model it is possible to define a 
'work term', contrary to what happens in the payment system models, as 
underlined in [20j. In order to show that the introduction of this term is 
consistent in this model, we define the monetary analogue of a Carnot heat 
engine, and show that the derived cycle has a consistent interpretation in 
terms of monetary policy. 

The monetary Carnot engine is ideal since isothermal transformations re- 
quire that the central bank is able to expand the monetary base without 
minimally affecting agents attitude regarding credit creation. For this rea- 
son, the monetary Carnot engine, as in thermodynamics, provides an upper 
bound for the effectiveness of monetary policy. This argument becomes ap- 
parent in the fractional reserve model where, taking into account (j43p . for 
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(59) 



isothermal transformations we have: 



1 \ F (I \V' 

v V 

V -V = (62) 
r r 

The r.h.s. of the last equation represents the variation of money supply, 
which is set to be equal to the variation of the monetary base. 
As underlined above, in adiabatic expansions an increase (decrease) of credit 
is associated with a decrease (increase) of the monetary base. In the case of 
real world monetary policy, this is a reasonable outcome. In fact, monetary 
easing is usually associated to a contraction of credit and vice versa, since 
the former is a reaction to the latter. 

In the monetary Carnot engine we obtain for L, with temperatures Tc, T^, > 
Tc, the following expression for L 

L = in-T,)iS'-S). (63) 

In thermodynamics the performance parameter ij of an engine working be- 
tween two sources with temperatures Tc,Tfi, > Tc is defined as 

\Qh\ 

where L is the work made by the engine and \Qh\ is the heat absorbed from 
the source at temperature T/j. For a Carnot engine we have 

l-|^<l-^ = r?, (65) 

\Qh\ J-h 

where the equality holds for quasi-statical transformations. The Carnot the- 
orem states that the performance of a given heat engine working between 
two temperatures Tfi,Tf. can never exceed the one given by a Carnot engine 
working with only reversible transformations. In our credit market model, 
T/j, Tc define two levels of average credit, with the policy performance factor 
7] given by 

where L is the work done by the central bank through monetary policy 
and \Ch\ is the overall credit variation. The policy performance factor can 
be interpreted as as measure of the effectiveness of monetary policy. By 
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maximizing ry, in other words, we maximize the impact of monetary pohcy 
over credit creation. As explained above, this impact is maximal if the 
central bank is able to follow the monetary Carnot cycle. Hence, instead of 
(j65|) we have: 

1_M<1_Z£ ^ \9A>Il, (67) 
\Ch\ Th \Ch\ Th 

For a constant number of agents we get: 

\Cc\ nic , s 

' ' > — (68) 



\Ch\ ruh 

Equation (j68p is an inequality relating measurable economic quantities. In 
particular, it states that the evolution of market attitudes, as given by the 
r.h.s. ratio, provides an upper bound for the effectiveness of monetary policy. 



7 Results and final remarks 

In this paper we built ab initio the statistical ensemble representation of two 
simple economic models. The starting point was the assumption of the con- 
servation of money, defined as equivalent to the monetary base, for a system 
composed by a given number of agents |1H \T7[ I13j . By using standard sta- 
tistical techniques, we wrote the economic analogue of all the main formulas 
of statistical mechanics where the role of the Hamiltonian was played by the 
money supply function M. In our context, the complete Hamiltonian for- 
malism of the ordinary statistical mechanics is not necessary. We only need 
the fact that there exists a function (the money supply function) which can 
be considered constant over a suitable time interval. In fact, the coordinates 
we use, i.e {x,y), are not conjugate. The analogy between the thermody- 
namics of a physical system and that of a credit market is based on a sound 
economic principle, i.e. the agents cannot 'print' but only exchange money. 
Only a central bank can print money. Once the identification between ther- 
modynamical quantities and economic ones is performed, we can take into 
account 'economic' transformations in which money can be exchanged be- 
tween different economic systems or work can be performed under the form 
of a credit expansion induced by monetary policy. 

Although real economic systems are typically not at equilibrium, this ap- 
proach can be used in principle whenever we can introduce a given economic 
conserved quantity. In particular, the study of equilibrium configurations 
can be useful to investigate the spontaneous evolution of those economic 
subsystems in which the conservation principle may be reasonably assumed 
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Table 1: Temperature T of the economic systems described by the money 
function M and by the ahowed credit and debt. 



as vahd. These are particular examples of markets in which, with some 
approximation, supply can be taken as fixed. As we have argued in the 
introduction, the credit market is an instance of this category inasmuch as 
the monetary base is determined exogenously by the central bank. 
In the context of our formalism we recovered the result of Dragulescu et 
al. [m [I2l [13] for a system described by the money function (fT5]) . In this 
case we obtain the same result for the economic temperature present in 
[llj . with and without the presence of debt (see formula (|27p . We stress 
that this result is obtained as a computation performed by employing the 
well-known mechanical statistic approach for the microcanonical ensemble, 
while the result of [TT] is obtained by fitting with an exponential Boltzmann 
distribution the result of a numerical simulation. Moreover we believe it 
is worth noticing the linear dependence on the variables Xj of the money 
function. This is the most striking difference from the thermodynamics of 
a physical system which involves a kinetic term with depends quadratically 
on Xi . We extended our analysis by considering generalization of the money 
function given by ([2]) by introducing the variables yi whose aim is to describe 
agents with different classes of assets and liabilities. We considered different 
models, without constraint on E ^« °^ Ey«> whose results are conveniently 
summarized in table ([1]) where we report M and T, which have different 
interpretations in the different models presented, and the presence of debt 
and/or credit. As a final application of our formalism, we explicitly show in 
the appendix which money function M is needed in order to recover a power 
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law Pareto distribution. For a discussion on this issue see also |21] . We leave 
for future work the application of our approach to other models from which 
different probability distributions are derived (see [22] for a review). 
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Appendix 

In [TU] it has been shown how if the following law holds 

dL 

for an agent i with income or wealth Ij, under the hypothesis that is 
randomly distributed around zero, then there exists a conserved quantity 

N 

Y = ^Hli). As a consequence, the Pareto power law [23] arises. Thus 

i=l 

we can use the formalism developed in this paper to study the statistical 
mechanics and thermodynamics which reproduces the Pareto law for the 
income. To this aim we used the canonical ensemble to describe the micro- 
scopic states of the system. The calculations are similar to the ones of the 
previous sections but now the conserved quantity is y. In this case instead 
of ([20]) we get 

-/^n.(^)^. 

where T is the economic temperature relative to the conserved quantity Y 
and Tmax is interpreted, as we see below, as the maximum for the tempera- 
ture allowed by the system. The integral ()70p is adimensional. As a result, 
we have a relation involving T, N and Y, i.e. the mean value of Y. As an 
example, we consider again a 'free' income agent function where Ij = Xj. 
Similarly to (I2ip . we have Z = e~T ^ F = Y — TS {S is the entropy func- 
tion) and {dF / dT)y ^ = —S. Furthermore, Pareto law is no more valid for 
small incomes li (for an empirical evidence see [10] and references therein) 
and as a consequence the integral in (j70p must be calculated with a lower 
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bound, i.e. Xj > J > 0. In this model we do not consider the variable yi in 
/. Finally, under the condition T < T^ax ensuring the convergence of the 
integral, we get: 



S = N ln{JVT) + N-N ln(r„,,. - T) + - 

J r 



NT 



(71) 



max 



-T 



AT ~ ' ^"■^ It I t — t' 

\ max / ^ max 



(72) 



For T << Tr 



max 1 



expression (j72]) becomes: 




max 



In 



( 



max 



J 



+ o{T). 



(73) 



Note that in this case, conversely to the analogous free case with M{xi), 
the temperature is not merely Y /N. Note that in general J < T^ax- As a 
consequence the second terms in the right hand side of ()73p gives a positive 
contribution to the temperature. Furthermore, as pointed in the study 
of quantities such as TdS/dT can give suitable information on possible phase 
transitions of an economic system. 

In our case, the calculation TdS/dT indicates an economic transition phase 
when T = Tmax- In fact for T = Tmax the expressions for Y ,S,F^Z diverge. 
This implies that the system can be in equilibrium only for T < Tmax- 
This simple example shows how our formalism can give some predictions 
concerning the possible macroeconomic changes of a given economic systems 
representing a country. 
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